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Abstract 

In LEP searches for the neutral CP-odd scalar A of a multi-Higgs doublet 
model, experimenters have searched for Z — > h°A° (where h° is the lightest CP- 
even scalar). No model- independent limit on the A mass can be deduced from 
present data if m^o > mz- In this paper, we compute the rates for Z — ► A G } A°uu 
and e + e _ — > A°A°Z. Assuming that no light CP-even neutral scalars exist, the 
branching ratio for Z — > A® A^vv is found to be less than 1.4 x 10 -8 . At higher 
e + e _ center-of-mass energies, er( e+e- -> A°A°Z) peaks at 0.5 fb. A comparison 
with other direct searches for the A is briefly considered. 
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1. Introduction 



The simplest extension of the scalar sector of the Standard Model is a model 
with two complex scalar Higgs doublets (with a CP-conserving Higgs sector) [1]. 
This model possesses five physical scalars: two neutral CP-even scalars (hP and H° 
with m/jo < mjjo), one neutral CP-odd scalar (A ) and a charged pair of scalars 
(i? ). In a generic (CP-conserving) two-Higgs-doublet model, the four masses m^o, 
mjjo, m^o and mjj± are independent parameters. Moreover, the scalar couplings 
to vector bosons and fermions depend on two additional independent parameters 
of the model: the ratio of vacuum expectation values (tan/9) and the CP-even 
Higgs mixing angle (a). This makes the experimental search for the various Higgs 
scalars a complicated task. 

The fact that the charged scalars have not been observed in Z decays puts 
an unambiguous lower bound on their masses, m#± > 41.7 GeV [2,3]. The Higgs 
search is more complicated in the case of the neutral scalars. The Z-boson does 
not couple to a pair of identical scalars but only to h°A° and H°A°. The relative 
coupling strengths are cos(/9 — a) and sin(/3 — a) respectively, which in the general 
two-doublet model is independent of the Higgs masses. Thus, if m^o > mz or if 
| cos(/9 — a) | <C 1 and m#o > mz, then associated production of A in Z decay 
could not have been observed even if A were massless! Experimental limits on A 
that do not depend on the other Higgs mass parameters of the theory are meager. 
There are no official limits in the 1992 compilation of the Particle Data Group [4] . 
Presumably, combining results from searches for K — > nA°, B — > KA°, ip — > A°7 
and T — > A ^ would rule out very light CP-odd scalars [5,6]. However, no general 
m^o bound above a few GeV exists today. 

Sharper bounds on m^o do exist in the minimal supersymmetric extension 
of the Standard Model (MSSM). The Higgs sector of the MSSM is a two-Higgs- 
doublet model in which the Higgs parameters are highly constrained [7]. As a 
result, one finds tree-level relations among the various Higgs sector parameters 
which reduce the six independent parameters mentioned above to two independent 
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parameters. These relations are modified when radiative corrections are incorpo- 
rated [8-11]. Nevertheless, the MSSM constraints permit a much more powerful 
set of experimental bounds to be obtained by the non-observation of Z — > hPvv 
and Z — > h°A°. The present LEP bound on A in the context of the MSSM is 
m^o > 20 GeV (assuming tan/? <: 1) [3,12]. 

In this paper, we consider the possibility of excluding a light A in the general 
(CP-conserving) non-supersymmetric two-Higgs-doublet model. We examine the 
extreme case in which the CP-odd neutral scalar A is light, but all other scalars 
are heavy (say, with masses above the e + e~ collider center-of-mass energy)* We 
then pose the following question: is it possible to find (or exclude) the A in a 
completely model-independent way? Since the other Higgs scalars are assumed 
to be heavy, the A must be produced either through its couplings to fermions 
or vector bosons. The Higgs-fermion interaction is model-dependent; the A°ff 
coupling strength depends on tan (3. On the other hand, the four-point couplings 
of A A to W + W~ and ZZ and their respective coupling strengths are fixed by the 
gauge invariance of the theory. Thus, we shall restrict our attention in this paper 
to A production processes that only involve the A couplings to vector bosons. 

In section 2, we study the decay mode Z — > A Q A Q vv shown in Fig. I. The vv 
pair arises from decay products of a virtual Z-boson. Our results can be easily 
extended to other pairs of light fermions arising from the virtual Z. If the heavier 
Higgs bosons of the model are sufficiently heavy, then the Feynman diagram of 
Fig. 1 is the only relevant tree-level diagram. If m^o = 0, we can compute the 
decay rate explicitly by analytically integrating over the four-body phase space. 
This provides an upper limit for the decay rate. We have also carried out the more 
general computation for m^o ^ where the phase space integration must be carried 
out numerically. In section 3, we compute the cross section for e + e~ — > A°A°Z. 
This process proceeds via a diagram similar to the one shown in fig. 1, except 
that the roles of the real Z and virtual Z* are reversed. The end results of both 

* The case where m^o <C raz and all other Higgs states are heavier than the Z is not possible 
in the MSSM. 
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calculations are disappointing; the respective branching ratio and cross-sections are 
too small to be detected at LEP (and LEP-200). Finally, in section 4, we compare 
our result with other direct searches for the A . 



2. Computation of V(Z -> A°A°vis) 

In this section, we compute the contribution of Fig. 1 to Z — > A Q A Q vv. The 
rate depends on the ZZAA coupling which in the two-Higgs doublet model is given 
by 



Note in particular that this coupling does not depend on the Higgs sector parame- 
ters a and (3. Thus, the decay rate we compute will depend only on one unknown 
parameter — the mass of the A . 

The amplitude for Z(p) — > A°(ki)A°(k2)i'(ks)^(k4 : ) (with the four-momenta 
indicated in parentheses) in unitary gauge is 



•ka ig2 nf ^(Sy-v- ^v/m 2 z ) _ ( -ig 

%M = {P) U{h) { 4^% ' 7 ^-kM**) , (2) 



where n = k% + k^ is the four-momentum of the virtual Z. The decay width is then 



T(Z - A°A°uu) = ^ f d " h d " h 1 

96(27r) 8 m| cos 6 9 W J 2E 3 2E A (2k 3 • h - m\f 

(3) 



x J [m 2 z (k 3 ■ h) + 2(k 3 ■ p)(h • P)} *< 4 >(p -h-k 2 -k 3 - k 4 ) , 



where we have included a factor of \ to account for identical A bosons in the final 
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state and a factor of i for the Z spin-average. To integrate over d^k\d^k 2 we use 



/ 



Sk x d 3 k 2 
2E\ 2E 2 



5 (4 \N-k 1 -k 2 )x< 



1 

kijj,kii 



_ 7T p 1 /^ 

— 2-"-4 ' 

= ^Rl /2 (R 4 N 2 g^ + 2R- 2 N fx N v ), 



(4) 



where TV is an arbitrary four-vector and 



R n = 1 — 



nm 



AO 



N 2 



(5) 



Applying this result to eq. (3), we find 



T(Z -> AAvv) 



d 3 ks d 3 k 4 



4:m 2 A0 



3(47r) 7 m| cos 6 % 7 2£ 3 2E 4 \ (p - k 3 - h) 2 
m 2 z (k 3 -k 4 ) + 2(k 3 ■ p){h • p) 



1/2 



X- 



[2(ks ■ k 4 ) - m z 



2 12 



(6) 



It is convenient to work in the Z rest frame: 

p = (m z ; 0,0,0); fc 3 = £ 3 (1;0,0, 1); k 4 = E 4 (l; sintf, 0, cos#), (7) 

where ^3 and E 4 are the neutrino energies. We introduce the following scaled 
kinematic variables 



4ml 
v = — w 



ID 



2 ' 



1 - cos 6 2E 3 2E 4 

— n — , y = — , z= — 



mz 



mz 



(8) 



f Here, we follow the techniques of ref. [13]. 



Then 



T(Z -> AAvv) 



g 6 m z 



l-v 



3-2 17 7T 5 COS 6 % 



dz 
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y dw + J dy J dw 

l-v-z (y+z-l+v)/yz 



x yVq + rg) ( , _ 
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y/ 2 

1 — y — z + yzw / 



(9) 



For a massless A (v = 0), the square- root factor in eq. (9) is unity, and the 
remaining integrals can be carried out explicitly. The result is 



T(Z -> AAvv) 



g 6 mz 



3-2 17 vr 5 cos 6 % 



7T 

T 



349' 
72~ 



(10) 



It is unfortunate that the factor in brackets is rather small (roughly 0.09), thereby 
reducing the partial width by more than an order of magnitude over a naive initial 
estimate. The branching ratio is most easily displayed by normalizing our result 
to F(Z -> vv) = g 2 m z / (96vr cos 2 9 W ). We find 



T(Z -> A A 



« 2 (f 



349 ' 
72 > 



T{Z -> i/P) 256tt 2 sin 4 % cos 4 % ' 



(11) 



where a ~ 1/128 is the fine-structure constant evaluated at the Z mass. Noting 
that the branching ratio for Z decay to three generations of neutrinos is 20.2%, we 
end up with 

BR(Z -> A°A°vv) = 1.4 x 10~ 8 . (12) 



For m^o ^ we have carried out the integrations in eq. (9) numerically. The 
results are presented in Fig. 2. The branching ratio decreases from the result of 
eq. (12) by about an order of magnitude for m^o ~ 15 GeV and by two orders of 
magnitude for m^o ~ 25 GeV. 
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3. e + e — > A°A°Z at a Higher Energy e + e Collider 

In light of the extremely small branching ratios for Z — > A A vis obtained 
above, we extend our considerations to future Higgs searches at higher energy e + e~ 
colliders (LEP-200 and the next e+e~ linear collider (NLC) with y/s > 300 GeV). 
Consider first the search for the lightest CP-even Higgs (h°) at LEP. If hP is not 
discovered at LEP-I (say, via Z — > h°vv), then the search will be continued at 
LEP-II via e+e- -> h°Z. Once the light CP -even Higgs state is found, one can, in 
principle, place unambiguous bounds on the existence of A if e + e~ — > h°A° is not 
observed. On the other hand, suppose that h° is still not discovered. Is it possible 
that a rather light A could still elude experimental detection? 

It might be argued that it is theoretically unlikely that A is light while the 
lightest CP-even scalar is quite heavy. This prejudice is common among advocates 
of low-energy supersymmetry models [14]. In contrast, light CP-odd and heavy CP- 
even scalars are natural in technicolor models [15]. For example, many technicolor 
models possess a spectrum of light CP-odd pseudo-Goldstone bosons ("pseudos"); 
those states that are electrically and color-neutral would behave like the A of 
the two-doublet Higgs model. In contrast, the CP-even scalars would be spin 
bound states of techni-fermions and would almost certainly be considerably heavier 
than the light pseudos. Thus, until a CP-even Higgs scalar is discovered, it will be 
important to study methods for discovering the A in processes that do not involve 
the CP-even scalars. 

In analogy with the search for the CP-even scalar at LEP-200, we consider 
the mechanism for A A production shown in fig. 1, where the role of the Z and 
Z* are interchanged. That is, we consider e + e~ — > A°A°Z via s-channel Z*- 
exchange. Here, we assume that all other Higgs scalars are much heavier than 
the e + e~ center-of-mass energy y^s so that other contributing tree-level diagrams 
mediated by the heavier CP-even scalars can be omitted. In this case, only the one 
diagram contributes. The squared matrix element for e + e~ — > A°A°Z, averaged 
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(summed) over initial (final) spins is given by: 

\kj\2 g 6 [l-4:sm 2 e w + 8sm 4 e w ] f 2 , 

|A1 1- = Mm%co S >O w (s-m%Y> ^ + 4?>1 ' " z P2 ' ^ ' < 13) 

where pi, p2 and k z are the e~, e + , and Z four- momenta respectively. In evaluating 
the cross-section for e + e~ — > A°A°Z, we make use of eq. (4) in the integration over 
the three-body phase space. This leaves one last integral over the three-momentum 
of the Z. At this point it is convenient to work in the e + e~ center-of-mass frame. In 
this frame, the energy of the Z is E z = (s+m z — S2)/2y / s where S2 is the invariant 
mass of the A A pair recoiling against the Z [16]. From 4m 2 40 < S2 < (y/s — mz) 2 , 
we easily obtain the limits of integration for the variable E z . Introducing the scaled 
variables: 

4mi m z E z 

v = — x = -=, y=-f, 14 

we end up with: 
a(e + e - A°A°Z) - 96s2 [1 " 4 Sin2 % + 8 Sin4 ° w] 



x 2 — 2y 



1/2 



3 • 2 12 vr 3 cos 6 e w m 2 z (s - m 2 z ) 2 

\{l+x 2 -v) 

x f dy(y 2 - X y/ 2 (y 2 + 2x 2 )(l- T 

X 

(15) 

In the above formula, we have already integrated over the angular direction of the 
Z, and we have included the factor of \ due to identical A bosons in the final 
state. Once again, we can perform the final integral analytically when m^o = 0. 
The result is 

, + - 4 o 4 o 7 x _ £ 6 s 2 [l - 4 sin 2 6 W + 8 sin 4 W ] 

1 J S^^TrScose^mKs-ml) 2 (lg) 

x [(l-x 4 )(l + 16:r 2 + x 4 ) + 72:r 4 lnx] . 

It is perhaps surprising that in the limit of large s, these cross-sections approach 
a constant value rather than behaving like 1/s as expected for annihilation cross- 
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sections. This behavior arises from the production of a longitudinal Z in the 
final state. Moreover, had we included the other tree-level contributions involving 
the exchange of the CP-even scalars (both h° and H°), we would find that the 
leading 0(1) behavior of the cross-section cancels exactly, leaving the expected 
0(1/ s) behavior. We have checked explicitly that this cancellation does occur. 
This means that eq. (15) and (16) can only be accurate for y/s smaller than the 
CP-even scalar masses; for larger values of y/s, these formulae must overestimate 
the true results. 

We now turn to the numerical implications of the above results. In fig. 3, 
we plot the cross-section for e + e~ — > A°A°Z as a function of the center-of-mass 
energy y/s for three choices of m^o. We also exhibit the dependence on m^o for 
possible energies of LEP-200 and the NLC (with y/s = 500 GeV). The cross- 
sections are disappointing. At LEP-200, optimistic projections suggest yearly data 
accumulations of about 500 pb _1 . This would imply less than one A°A°Z event 
per year, even before detection efficiency and Z branching ratio factors are taken 
into account. At the NLC, typical design luminosities yield between 20 and 50 fb _1 
of data per year [17]. In this case, at most a few A°A°Z events per year could 
be detected (after experimental triggers and cuts). Of course, such a search would 
only be critical if no h° were discovered at the NLC, as discussed at the beginning 
of this section. 

Finally, for completeness, we note that other production processes for A A 
should be considered at the NLC once the center-of-mass is sufficiently large. For 
example, one can consider e+e~ -> A°A°W+W- and e+e~ -> A°A°ZZ where 
the A A is emitted from one of the outgoing vector boson lines. In addition, 
for y/s > 500 GeV, vector boson fusion processes begin to become important. 
Thus, one should examine e + e~ — > A® A^vv (via W + W~ fusion) and e + e~ — > 
A°A°e+e- (via ZZ fusion). The cross-sections for all these processes are formally 
of higher order [(9(g 8 )] as compared to e + e~ — > A°A°Z considered above. However, 
the cross-section for the fusion processes increase logarithmically with y/s and 
eventually dominate the annihilation processes. The W + W~ fusion cross-section 
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is the largest and was computed in ref. 18. Unfortunately, even at ^/s = 2 TeV, its 
cross-section appears to be too small for detection at a future e + e~ linear collider. 

4. Discussion and Conclusions 

The branching ratio for Z — > A Q A Q h>v is probably too small to be observed 
even at a high luminosity Z factory. Including other possible decays of the virtual 
Z to fermion pairs can only enhance the branching ratio by a factor of 5. Even 
with 100% detection efficiency, it seems unlikely that the A can be detected. 
Our calculation has assumed the most general (CP-conserving) two-Higgs-doublet 
model with a light A and all other Higgs scalars very heavy. If the other scalars 
are not all heavy, there are additional diagrams that contribute to Z — > A® A^vv. 
We expect an enhancement of the Z — > A°A°i>i> decay rate to be at least an order 
of magnitude beyond the result of eq. (12), since there is no reason for a similar 
suppression by the factor 7T 2 /2 — 349/72 as found in eq. (10). After adding up 
all diagrams involving the other Higgs states of the model, we would expect a 
branching ratio perhaps as large as 10~ 6 . A more complete calculation (in the 
context of the MSSM) is currently in progress [19]. Such a decay rate may be 
observable at a future high intensity Z factory. However, a favorable branching 
ratio would exist only in limited regions of the parameter space. Going to a higher 
energy e + e~ collider does not improve the picture substantially. In this case, one 
would search for e + e~ — > A°A°Z. Based on our computations, we would conclude 
that the production cross sections are not large enough to permit experimental 
detection at LEP-200. At an even higher energy NLC with a data collection greater 
than 10 fb" 1 per year, at most a few signal events per year could be detected. Thus, 
an unambiguous model-independent limit on the A mass seems to be unattainable. 

It is instructive to compare the results of this paper with other proposals to 
search for A in rare Z decays. Two comparable decay processes have been studied 
in the literature: Z -> A A A and Z -> bbA°. L.-F. Li has estimated [20] the 
rate for Z — > A A A via a one-loop box diagram with a t-quark running around 
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the loop. He estimates BR(Z — > A A A ) ~ 10~ 5 , although his estimate is rather 
crude. If this estimate is not overly optimistic, then it clearly presents an excellent 
chance to obtain an experimental bound on m^o which is independent of the other 
Higgs masses of the model. However, there is still some model-dependence in this 
case, since the amplitude for Z — > A A A (via the t-quark loop) depends on the 
A°ti coupling. This coupling depends on tan (3 and could be suppressed. Of course, 
if h° is not too heavy, then the dominant contribution to Z — > A A A would arise 
from Z — > A°h°* in which the virtual h°* decays to A A . 

The decay rate for Z — > bbA° also depends on the model-dependent Higgs- 
fermion coupling. Djouadi, Zerwas and Zunft find BR(Z — > bb~A°) ~ 3 x 1CP 4 for 
rriAo = 10 GeV and tan/5 = 20 [21]. (This branching ratio scales as tan 2 (3.) Once 
again, we have the potential for obtaining limits on A independent of the masses 
of the other Higgs scalars of the theory. As before, the mo del- dependence of the 
Higgs-fermion coupling still remains. However, if we subscribe to the Model-II 
Higgs-fermion couplings of ref. [1] (this coupling pattern is the one that appears in 
the MSSM), then the A coupling to bb (ft) is proportional to tan/5 (cot j3). Thus, 
the latter two processes are complimentary. A combined search for Z — > A ^ ^ 
and Z -> bbA° at a future high intensity Z factory could then rule out the possible 
existence of a light A . At a higher energy e + e~ collider (such as the NLC), one 
can search for both e + e~ — > bbA° and e + e~ — > tiA° with similar effect. On the 
other hand, in the Model-I Higgs-fermion coupling of ref. [1], it is possible that 
both A couplings to bb and ti are suppressed. This leaves open the possibility 
that a light A could escape experimental detection for the foreseeable future. 
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FIGURE CAPTIONS 



1) Feynman diagram for Z — > A® A^vv. In the limit where all other Higgs scalars 
are much heavier than mz, other contributing tree-level diagrams (mediated 
by virtual heavy Higgs exchange) can be neglected. 

2) Phase space suppression of Z — > A Q A Q vv as a function of m^o. Tq is the 
corresponding decay rate for m^o = 0; at this point, BR(Z -> A°A°z/z/) = 
1.4 x 10~ 8 . 

3) Cross-section for e + e~ — > A°A°Z. We exhibit curves in two cases: (a) as a 
function of y/s for m^o = 0, 20 and 40 GeV and (b) as a function of m^o for 
v/i = 150, 175, 200 and 500 GeV. 
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